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Abstract 

The previously developed technique for evaluation of charge-transfer and electron-excitation processes in low- 
energy heavy-ion collisions [I.I. Tupitsyn et al, Phys. Rev. A 82, 042701(2010)] is extended to collisions of ions 
with neutral atoms. The method employs the active electron approximation, in which only the active electron 
participates in the charge transfer and excitation processes while the passive electrons provide the screening DFT 
potential. The time-dependent Dirac wave function of the active electron is represented as a linear combination of 
atomic-like Dirac-Fock-Sturm orbitals, localized at the ions (atoms). The screening DFT potential is calculated 
using the overlapping densities of each ions (atoms), derived from the atomic orbitals of the passive electrons. The 
atomic orbitals are generated by solving numerically the one-center Dirac-Fock and Dirac-Fock-Sturm equations 
by means of a finite-difference approach with the potential taken as the sum of the exact reference ion (atom) 
Dirac-Fock potential and of the Coulomb potential from the other ion within the monopole approximation. The 
method developed is used to calculate the K-K charge transfer and K-vacancy production probabilties for the 
Ne(ls 2 2s 2 2p 6 ) - F 8+ (ls) collisions at the F 8 +(ls) projectile energies 130 keV/u and 230 keV/u. The obtained 
results are compared with experimental data and other theoretical calculations. The K-K charge transfer and 
K-vacancy production probabilities are also calculated for the Xe - Xe 53+ (ls) collision. 

PACS numbers: 34.10,+x, 34.50.-s, 34.70.+e 
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I. INTRODUCTION 



Collisions of highly charged ions provide a unique tool for tests of relativistic and quantum electro- 
dynamics (QED) effects in the scattering theory yj-yl]. Investigations of such processes can also give 
an access to QED in supercritical fields, provided the total charge of the colliding nuclei is larger than 
the critical one, Z c = 173 (see, e.g., Ref. [4] and references therein). One of the most attractive way 
for indirect observation of the supercritical field created in the collision is to investigate the dynamics 
of inner-shell electrons that can be rather sensitive to the collision parameters. The most favorable con- 
ditions for studying the electron dynamics in the supercritical field regime correspond to the projectile 
energy of about the Coulomb barrier [5]. In case of U— U collision this requires the energy of about 
6 MeV/u that means a low-energy collision. One of the key processes in such collisions is the charge 
transfer of electrons (see, e.g., Ref. [6] and references therein). A systematic approach to relativistic cal- 
culations of the charge-transfer and electron-excitation probabilities in low-energy heavy-ion collisions 
was developed in our previous paper [7], where the consideration was restricted to collisions of H-like 
ions with bare nuclei. Since the experimental study of such collisions for high-Z systems is presently 
rather problematic, an extention of the method to collisions of highly charged ions with neutral atoms, 
that can be studied in experiments with the current GSI and future FAIR facilities fl8|,|9D, is needed. In 
this paper we present the desired extention and perform calculations for low-energy ion-atom collisions. 
To examine the approach we calculate the K-K charge transfer and K- vacancy production probabilities 
for low-energy collision of H-like F (F 8 +) and neutral Ne, the process which has been investigated both 
experimentally | j 1 Qfl and theoretically flliHlSQ. The calculations are performed at the F 8+ projectile 
energies 130 keV/u and 230 keV/u. We also evaluate the probabilities of the K-K charge transfer and 
K- vacancy production in the Xe - Xe 53+ (ls) collision at the projectile energy of 3.6 MeV/u. The latter 
processes are planned to be studied in the nearest future in experiments at GSI [8, ^]. 



The paper is organized as follows. In the section UTAI we describe the time-dependent one-electron 
equation in so-called active electron approximation [14] and the method for constructing the local 
Kohn-Sham potential induced by the passive electrons. The wave function of the active electron is 
expanded in terms of the Dirac-Fock and Dirac-Fock-Sturm basis functions, which are central-field 
4-component Dirac bispinors centered at the ions. The two-center relativistic Kohn-Sham equation in 
the finite basis set is briefly discussed in the section ITTBl The basis functions are obtained by solving 
numerically the atomic Dirac-Fock and Dirac-Fock-Sturm equations in an external field with a special 
choice of the weight function, as it was proposed in Refs. [ 1^1 . The external potential is spherically 



symmetric Coulomb-Hartree potential of the other ion (atom) taken in the monopole approximation. 
The basis set constructed in this way and the related calculation procedures are described in section HTCl 
Basic formulas for the K- vacancy production probability are given in section III Dl In section III El we 
present the results of the relativistic calculations of the Ne K-shell-vacancy production and K-K charge 
transfer probabilities in the Ne - F 8+ collisions as a function of the impact parameter b at the projectile 
energies 130 keV/u and 230 keV/u. In this section we also present the results of the neutral Xe K-shell- 
vacancy production and K-K charge transfer probabilities in the Xe - Xe 53+ (ls) collision. 



II. THEORY 



A. Dirac-Kohn-Sham equation in the active electron approximation 



In this paper we use so-called active electron approximation |14J] to describe the ion-atom collision. 
In this approximation, we consider only the active electron which participates in the charge transfer and 
excitation processes, while the other passive electrons provide a screening potential. In our calculations 
the screening potential is defined by the density functional theory (DFT) in the local density approxima- 
tion (LDA). In this approach the time-dependent wavefunction ip(r, t) of the active electron ^(r, t) is 
the solution of the relativistic time-dependent Kohn-Sham equation. In atomic units (h — m — e — 1), 
this equation is given by 

J^-fco^r,*). (1) 
Here is the two-center Dirac-Kohn-Sham Hamiltonian defined by 

h D = c(cx-p) + (3c 2 + V AB (r) , V AB (r) = V H [p\ + V xc [p] , (2) 

where c is the speed of light and a, (3 are the Dirac matrices. V B [p] and V xc [p] are the Hartree and 
exchange-correlation potentials, respectively. Both of them are the functionals of the electron density 
p(r). The Hartree potential V B [p] includes the electron-nucleus interaction and the electron-electron 
Coulomb repulsion Vc[p\: 

V H (r) = V* cl (r A ) + V* cl (r B ) + V c [p] , r A = r-R A , r B = r - R B , (3) 



where 



j3 j /w(r-') T . { , f , 3 , p(r') 



Kuci(r) = / dV , Vcip] = I d A r' , (4) 
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Pnuci(^) and p{r) are the nuclear and electron densities, respectively. The exchange-correlation poten 



tial V xc [p] was taken in the Perdew-Zunger parametrization Ill7h including the self-interaction correction 
(SIC). 

The electron density p(r), obtained with the many-electron wave function which is represented by 
a Slater determinant, is invariant with respect to the rotations in the occupied orbitals space. For this 
reason, to obtain the electron density we can use atomic-like orbitals, localized on the both centers (A 
and B). In this case the electron density, constructed from the orthogonal localized orbitals, can be 
represented as a sum of densities pa(^) and p B (r) which are localized on the centers A and B. This is 
not the case, however, if the orbitals localized on the different centers overlap and are non-orthogonal. 
The electron density, derived from non-orthogonal orbitals, is given by 

p(r) = J24 Pr (r) (S- 1 )..^ (r), (5) 

where ip^ i r ) are the atomic-like wavefunctions of the passive electrons and matrix S is the overlapping 
matrix. Note that the electron density is normalized on the number of the passive electrons, 

J d 3 r p{r) =N-1, (6) 

where N is the total number of electrons. The electron density p(r) can be divided into three parts, 

p(r) = p A (r) + p B (r)+p^ P \r), (7) 

if we split the summation over indices i,j into the sum over i,j e A, the sum over i,j E B and the 
remaining overlapping part. We can also split the overlapping density into two parts dividing the space 
into two regions (A) and (B), 



This can be done by the plane passing through the middle of the internuclear distance (see Ref. [7D for 
details). 

For simplicity, let us consider spherically-average values of the electron densities in each region, 
p A {r A ) = J dtt A p A {r) + p^ ovlp) (r) , p B (r B ) = J dVt 
This procedure does not change the normalization of the total electron density. 



(9) 



d A r p(r) = / d 6 r [p A (r A ) + p B (r B )] — N — 1. (10) 



As a result, the potential Vab{t) can be approximated by the sum of the spherically symmetric poten- 
tials of the two different centers, 

V AB {r) ~ V A [p A ](r A ) + V B [p B ](r B ) . (11) 

The overlapping densities must be taken into account especially for the short internuclear distances. 
Otherwise, the Pauli principle is violated and, as a result, the number of electrons on the Is shell of the 
united system can exceed 2. 

It should also be noted that the time-dependent wave function ip(r, t) of the active electron is or- 
thogonalized to the wave functions ij)f\r) of the passive electrons. This means that the transitions of 
the active electron to the states occupied by the passive electrons are forbidden in accordance with the 
Pauli principle. 



B. Two-center Dirac-Kohn-Sham equation 

The two-center expansion of the time-dependent wave function ip(r, t) can be written in the form 

il>{r,t)= ^C^t) <p a , a (r - Rait)) , (12) 

a=A,B a 

where index a = A, B labels the centers, index a enumerates basis functions at the given center, and 
^a,a{f — R a ) is the central-field bispinor centered at the point a. In what follows, the shorthand 
notations \j) = \(pj) = \(p a , a ) f° r states j = a, a are used. The expansion coefficients C aa (t) of 
the time-dependent wave function ip(r,t) can be obtained by solving the linear system of first-order 
differential equations 

k k 

where indices j and k enumerate the basis functions of both centers, and the matrix elements of H and 
S are 

H jk = (j \h D \k), S jk = (j | k) . (14) 
The matrix elements of T are given by 

T j k = i{J\^\k)=Tt j + i^S jh . (15) 
Obviously the matrix T is non-Hermitian, if the overlapping matrix S depends on time. 
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The functions ip a>a depend on time due to two reasons. First, the basis functions centered at the target 
and projectile nuclei move together with the nuclei. Second, the basis functions depend parametrically 
on the distance between the nuclei, since their radial parts are obtained from the radial equations, where 
for each center the potential of the other ion (atom) is included in the so-called monopole approximation 
(see section HTCl ). 

Calculations of the matrix elements Hjk, Sjk, and Tjk were considered in detail in Ref. QVU . 
C. Basis functions 

In our approach the basis set contains Dirac-Fock and Dirac-Fock-Sturm orbitals. The Dirac-Fock- 
Sturm orbitals can be considered as pseudo-states, which should be included in the basis to take into 
account the contribution of the positive- and negative-energy Dirac continuum. Both types of basis 
functions Lp aa are central field Dirac bispinors centered at the position R a (a = A, B) of the corre- 
sponding ion, 



r 



(16) 



r 

where P nK (r) and Q nK (r) are the large and small radial components, respectively, and k = 
(_iy+j+V 2 (j -f 1/2) is the relativistic angular quantum number. The large and small radial com- 
ponents are obtained by solving numerically the Dirac-Fock and Dirac-Fock-Sturm equations in the 
central field approximation. The radial Dirac-Fock equation is 

(/£ F + Kxt(r)) F anK (r) = e anK F anK (r) , F anK (r) = I PanK ^ I , (17) 

where h® F is the radial Dirac-Fock Hamiltonian of ion a (a = A, B), F anK {r) is the two-component 



radial wave function, and V ext (r) is a local external potential. The explicit form of the radial 



Dirac-Fock 



equation and the description of the corresponding computer code are presented in Ref. H18H . The radial 
components of the Dirac-Fock-Sturm orbitals Tp nKm , which we denote by F nK (r), are the solutions of 
the generalized Dirac-Fock-Sturm eigenvalue problem, 

(h% F + V ext (r) - e anoK ) F anK (r) = X anK W aK (r) F anK (r) . (18) 

Here \ anK can be considered as the eigenvalue of the Dirac-Fock-Sturm operator and W aK (r) is a 
constant sign weight function. The energy e anoK is fixed in the Dirac-Fock-Sturm equation. If the 
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weight function W(r) — > at r — > oo, all Sturmian functions have the same asymptotic behavior 
at r — > oo. It is clear that for X anK = the Sturmian function Tp nKm coincides with the reference 
Dirac-Fock orbital (p a n K- hi our calculations we use the following weight function 

In contrast to 1/r, this weight function is regular at the origin. It is well-known that the Sturmian oper- 
ator is Hermitian. It does not contain continuum spectra, in contrast to the Dirac operator. Therefore, 
the set of the Sturmian eigenfunctions forms a discrete and complete basis set of one-electron wave 
functions. 

The external central-field potential V ext (r) in equations (fT7l) and (fT8l) is arbitrary, and, therefore, it 
can be chosen to provide most appropriate Dirac-Fock and Dirac-Fock-Sturm basis orbitals. At small 
internuclear distances the wave function of the atomic electron experiences also the strong Coulomb 
field of the other ion (atom). To the leading order this effect can be taken into account by including 
the Coulomb-Hartree potential of the second ion (atom) as the external potential V cxt (r) within the 
so-called monopole approximation. For instance, the external central-field potential V^ t (r) is given by 

KxtW = VZn(r) = ^ J dn A V*{r - R AB ) . (20) 

where V^ on (r) is the spherically- symmetric part of the reexpansion of the Coulomb-Hartree potential 
V§ ( r — Rab) °f the ion B with respect to the center A and Rab is the internuclear distance vector. 

Calculations of the two-center integrals with the basis functions employed were described in detail 
in our previous paper [7]. The time-dependent Dirac-Kohn-Sham equation for the active electron is 
solved using the two-center basis set expansion. The expansion coefficients are determined employing 
the direct evolution (exponential) operator method |7|], which is more stable compared to the others, 
such as, e.g., the Crank-Nicholsen propagation scheme [19] and the split-operator method [20]. To 
obtain the matrix representation of the exponential operator in the finite basis set one has to diagonalize 
the generalized complex Hamiltonian matrix at each time step. Since our basis set is not too large, the 
diagonalization procedure is not too time consuming. 



D. Charge-transfer and vacancy-production probabilities 



The amplitudes of the charge transfer and excitations to different bound states of the projectile and 
target ions are calculated by projecting the time-dependent wave function of the active electron onto the 
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atomic Dirac-Fock orbitals of the projectile and target. The corresponding calculations for collisions of 
H-like ions with bare nuclei were described in detail in our work [7]. That is why here we restrict our 
consideration only to the new features of the calculation procedure that occur for atom-ion collisions 
within the one active electron approximation. 

Consider the collision of a neutral atom A (target) with a hydrogenlike ion B (projectile). We 
assume that before the collision the active electron occupies the Is state of the target with spin up (in 
the relativistic case, with the total angular moment projection jx = 1/2) and the passive electron of 
the hydrogenlike ion occupies the Is state with spin down. In what follows, we are interested in two 
processes: the K-K charge transfer and the K-shell vacancy production. 

Let (Is) is the probability to find the active electron in the Is state of the target after the collision, 
Pb(1s) denotes the probability to find one active electron in the Is state of the projectile or, in other 
words, the probability Pk-k of the K-K shell charge transfer of one electron. To obtain the probabilty 
P vac of the K-shell vacancy production, we introduce the probabilities P(Ei) and P{E 2 ) of the events 
Ei and E 2 , when a hole is created in the Is state of the target with spin up and down, respectively. The 
probabilities of these events are defined by 

P{E X ) = 1-Pa(1s) 

(21) 

P{E 2 ) = 1 - P A (ls) - P B (ls) 

Assuming the events E\ and E 2 are independent, the probability of production at least one hole in the 
Is state of the target is given by 

P vac = P{Ei) + P(E 2 ) - P(Ei)P(E 2 ) = 1 - P A (ls)(P A (ls) + P B (ls)) . (22) 

We note that, since the sum of the probabilities P^(ls) and Pg(ls) is less than 1, the vacancy production 
probability satisfies the condition < P vac < 1. It should also be noted that the K-shell vacancy 
production defined by Eq. (1221) includes the production of two holes in the target K shell. 



E. Results of the calculations and discussion 



To test the approach we have studied the Ne - F 8+ (ls) collision for low energies where experimental 
and nonrelativistic theoretical results are available [| 1 QM 1 311 . In this case the nuclear charge numbers are 
rather small and, therefore, relativistic effects are negligible. We stress, however, that our approach can 
be directly applied to heavier systems where the relativistic effects become stronger or even dominant. 



In Fig. Q] we present the results of our calculations for the probabilities P(b) of the Ne K-shell- 
vacancy production (solid line) and of the K-K-shell charge transfer (dotted line) as functions of the 
impact parameter b for the Ne-F 8+ (ls) collision at the projectile energy of 230 keV/u. The results for 
the Ne K- vacancy production are compared with experimental values (circles) IllOn and with theoretical 
results obtained by Fritsch and Lin (dashed line) [11] and by Thies et al. (dash-dotted line) 1130. It can 



be seen that our results are in perfect agreement with the experimental ones. 
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FIG. 1: The results of the present calculations for the probabilities P(b) of the Ne K-shell- vacancy production 
(solid line) and of the K-K-shell charge transfer (dotted line) as functions of the impact parameter b for the Ne- 
F 8+ (ls) collision at the projectile energy of 230 keV/u. The circles indicate experimental results by Hagmann 



et al. 



10]. The dashed and dash-dotted lines present theoretical results by Fritsch and Lin H 1 ill and by Thies et 



al. II 1311 . respectively. 



The related results for the Ne-F 8+ (Is) collision at the projectile energy of 130 keV/u are presented in 
Fig.|2l where the same notations as in Fig.[T]are used. We note that our theoretical results are in a good 
agreement with the experimental ones at small impact parameters. However, in contrast to Fig. I, the 
agreement is not so good for medium and large impact parameters, although the theoretical predictions 
for the maximum and minimum positions agree rather well with the experimental ones. 

As one can see from Figs. [Hand (2] for both energies at large impact parameters the K- vacancy 
production is mainly determined by the K-K-shell charge transfer, which is indicated by the dotted line. 
The difference between the K-vacancy production and the K-K-shell transfer probabilities at small 
impact parameters is due to the contribution from the charge-transfer excitation into the 2s, 2p, and 
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higher vacant states of the projectile. This is in accordance with the experimental results of Ref. Bill , 
where the K-vacancy production in the Ne - F 6+ ((ls) 2 2s) collision was studied. The calculation of the 
latter process is currently under way and will be published elsewhere. 

1 r 
P(b) 




0,6 0,7 0,8 

b [a.u.] 



FIG. 2: The results of the present calculations for the probabilities P(b) of the Ne K-shell-vacancy production 
(solid line) and of the K-K-shell charge transfer (dotted line) as functions of the impact parameter b for the Ne- 
F 8+ fls) collision at the projectile energy of 130 keV/u. The circles indicate experimental results by Hagmann et 



al. H 1QTI . The dashed and dash-dotted lines present theoretical results by Lin et al. (taken from Ref. [ 10]) and by 



Thies et al. H 1 3TI . respectively. 



In this work we also performed the related calculations for the Xe - Xe 53+ (ls) collision at the 
projectile energy of 3.6 MeV/u. The experimental study of this process is planned at GSI (Darmstadt) 
[8i|9|]. The probabilities P(b) of the Xe K-shell-vacancy production and of the K-K-shell charge transfer 
as functions of the impact parameter b are plotted in Fig. [3j The solid and dotted lines represent the 
vacancy production and the charge transfer, respectively. For comparison, in the same figure we display 
the K-shell-vacancy production for the Xe 53+ (ls) - Xe 54+ collision that is indicated by the dashed 
line. We note again that at large impact parameters the K-vacancy production is almost completely 
determined by the K-K-shell charge transfer. It can be also seen that in the case under consideration the 
screening effect is rather small. 

To investigate the role of the relativistic effects we performed the same calculations for the Xe- 
Xe 53+ (ls) collision in the nonrelativistic limit by multiplying the standard value of the speed of light 
by the factor 1000. The obtained relativistic and nonrelativistic results are presented in Fig. [4] As one 
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FIG. 3: The probabilities P(b) of the Xe K-shell-vacancy production (solid line) and of the K-K-shell charge 
transfer (dotted line) in the Xe-Xe 53+ (ls) collision as functions of the impact parameter b. The dashed line 
indicates the K-shell-vacancy production for the Xe 53+ (ls)-Xe 54+ collision. 

can see from the figure, the oscillatory behavior of both curves is the same but the nonrelativistic curve 
is shifted toward higher impact parameters. 



III. CONCLUSION 



In this paper the method that was previously developed for evaluation of the electron-excitation and 
charge-transfer processes in collisions of H-like ions with bare nuclei has been extended to collisions 
of ions with neutral atoms. The extention is based on the active electron approximation, in which 
the interaction of the active electron with the passive electrons is accounted for by the screening DFT 
potential. 

The method developed has been applied to evaluate the K-vacancy production and the K-K charge 
transfer in the low-energy Ne - F 8+ collision. The results of the calculation are compared with available 
experimental data and with theoretical calculations by other authors. The influence of the relativistic 
effects on the K-vacancy production probability is investigated for the Xe - Xe 53+ collision. It is 
demonstrated that the relativistic and nonrelativistic probabilities as functions of the impact parameter 
exhibit the same oscillatory behavior at low energies but the relativistic curves are shifted toward lower 
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Xe-Xe^ + (ls) 3.6 MeV/u 




b [a.u.] 

FIG. 4: The probability P(6) of the Xe K-shell-vacancy production in the Xe-Xe 53+ (ls) collision as a function 
of the impact parameter b. The solid and dashed lines present relativistic and nonrelativistic results, respectively. 

impact parameters compared to the nonrelativistic ones. 

In our further investigation we plan to continue calculations of low-energy heavy-ion collisions that 
are of interest for current and nearest future experiments at GSI and FAIR facilities in Darmstadt. Spe- 
cial attention will be paid to the critical regime, when the ground-state level of the united quasimolecule 
dives into the negative-energy Dirac continuum. 
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